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We present here the mass excesses, binding energies, one- and two- neutron, one- 
and two- proton and Q-particle separation energies of 6727 nuclei in the ranges 
4 < Z < 120 and 8 < A < 303 calculated in the infinite nuclear matter model. 
Compared to our predictions of 1999 mass table, the present ones are obtained using 
larger data base of 2003 mass table of Wapstra and Audi and resorting to higher 
accuracy in the solutions of the ry-differential equations of the INM model. The local 
energy ry's supposed to carry signature of the characteristic properties of nuclei are 
found to possess the predictive capability. In fact r/-systematics reveal new magic 
numbers in the drip-line regions giving rise to new islands of stability supported 
by relativistic mean field theoretic calculations. This is a manifestation of a new 
phenomenon where shell-effect overcomes the instability due to repulsive components 
of the nucleon-nucleon force broadening the stability peninsula. The two-neutron 
separation energy-systematics derived from the present mass predictions reveal a 
general new feature for the existence of islands of inversion in the exotic neutron- 
rich regions of nuclear landscape, apart from supporting the presently known islands 
around ^^Na and ®^Ti. The five global parameters representing the properties of 
infinite nuclear matter, the surface, the Coulomb and the pairing terms are retained 
as per our 1999 mass table. The root-mean-square deviation of the present mass-fit 
to 2198 known masses is 342 keV, while the mean deviation is 1.3 keV, reminiscent 
of no left-over systematic effects. This is a substantive improvement over our 1999 
mass table having rms deviation of 401 keV and mean deviation of 9 keV for 1884 
data nuclei. 



* This is a briefly enhanced version of the article just published in Atom. Data and Nucl. 
Data Tables 98(2012)213-719. 
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I. INTRODUCTION 

Mass formulas always occupy the center-stage in the research arena of nuclear physics. 
Traditionally mass formulas have been developed using the two main features of nuclear 
dynamics, namely, the liquid drop and the shell features. The liquid drop has been the 
mainstay of nuclear physics since its introduction by Weizsacher and Bethe in 1930s. The 
liquid - supposed to represent the nuclear matter composed of the nucleus - is a classical one, 
although in reality it is a quantum liquid consisting of interacting many-fermionic system. 
At macroscopic level the classical liquid picture has been most useful as a reference model. 

nn 

However it has given rise to some discomfitures, the notable being the tq — paradox\]l |2| 
i.e., the discrepancy between the density determined by it through the fit to nuclear masses 
and the actual density measured through electron scattering on heavy nuclei. In the infinite 
nuclear matter (INM) modeljs-G], the classical liquid drop is replaced by an INM sphere 
characterizing the interacting many-fermionic liquid. Further the model is based on the 
Extended Hugenholtz-Van Hove theorem Q] of many-body theory, and therefore takes 
into account the shell and the liquid-drop features nonperturbatively. 

The first mass table based on this model was published jsl in 1988 with predictions of 
3481 nuclei. The far-off nuclei in the drip-line regions could not be included because of 
accumulated error resulting from imprecise solution of r^-equation. Further development in 
the INM model was achieved by using better definition of Fermi energies leading to perfect 
decoupling 8^, M of the finite-size terms like surface and Coulomb from the infinite-nuclear- 

nn 

matter terms characterized by volume and asymmetry. This led to the determination [8|,l9j of 
the saturation properties of nuclear matter i.e., the density, the bindingenergy per nucleon 
and the incompressibility from nuclear masses and thereby resolvingj^, ^ the tq — paradox. 
Thus the basic objective of the INM model was achieved. Further for the successful long- 
range prediction of nuclear masses, an interactive network extending over the entire nuclear 
chart was devisedj^ for better determination of the local energy 77. This procedure yielded 
multiple values of rj for a given nucleus whose ensemble- average was carried out to obtain 
its unique value. With these improvements, reliable prediction of nuclear masses up to 
drip-line could be made and the second mass table containing the data for 7208 nuclei was 
publishedje] in 1999. 

Up to 1999 our concern was directed towards the long-standing problem of the determi- 
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nation of the saturation properties of nuclear matter and its incompressibility from nuclear 
massesjs, Once this objective was fulfilled we concentrated on the study of properties 
of local energy rj, a crucial element in the INM model, t] embodies all the characteristic 
properties of a nucleus, and therefore carries its finger-print. The r^'s determined in our 
1999 mass predictions were analyzed as a function of neutron number N for a given proton 
number Z revealing lOj Gaussian structure for the well-known closed-shell nuclei in the val- 



ley of stability and similar pattern for some nuclei in the far-off drip-line regions [ll|]. On 
the basis of such revelations, new neutron magic numbers 100, 150, 164 and proton magic 
number 78 and weakly proton magic number 62 and 90 were predicted HI] designating new 
islands of stability around ^^"^Sm, "^"^^Pt and "^^^Th reaffirmed through microscopic study H] 
in the Relativistic Mean Field (RMF) theory and Strutinsky Shell-Correction calculations. 
This is suggestive of a new phenomenon where shell stabilizes the instability due to repulsive 
components (triplet-triplet, singlet-singlet) of nucleon-nucleon interaction, complementary 
to the phenomenon of fission isomers and super-heavy elements in which repulsive coulomb 
instability is overcome by the same. While the latter elongates the stability peninsula, the 



former broadens it. In another development, comparative study 12 



14| of the predictive 



ability of different mass formulas carried out with respect to a reference mass formula shows 
random divergences without any common trend. This raises the question about predictive 
ability of mass formulas and their possible use as guidelines for theoretical and experimental 
study in unknown regions. Taking the masses calculated in the microscopic study employ- 
ing RMF theory as reference masses, it is shown 15|| that the divergence disappears yiel ding 
common trend for the prediction of all mass formulas . Similar trend is also observed |l5| 
when masses predicted in our 1999 mass table are used as reference. This success could 
be attributed to the very structure of the INM model equations. It is identified that this 
mass formula is written down in terms of differential equations. It is a well recognized fact 
that predictive ability is intrinsic in any theory formulated in terms of differential equations 
like that of Newton, Schrodinger and Dirac etc. involving rate of change of physical vari- 
ables. This feature of the INM model was not specifically recognized or perceived so far. 
We had implicitly articulated that r/'s of unknown regions were obtained by some kind of 
extrapolation, which was somewhat misleading. 

In view of the above developments after 1999, and quite importantly availability of large 
number of new measurement of mass of unknown nuclei in the intervening period, prompted 
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us to make a fresh and probably the final mass table. Addition of 314 newly measured data 
to the already used 1884 masses in our 1999 mass table, and more precision solution of rj- 
differential equation being possible has resulted in the reliable prediction of masses covering 
up to drip-line regions with rms deviation of 342 keV and mean deviation of 1.3 keV. This 
is a substantive improvement over our earlier result of 1999 for 1884 mass data with rms 
deviation of 401 keV and mean deviation of 9 keV. 



II. THE INFINITE NUCLEAR MATTER MODEL 

Although the INM model is well-known and described in several papers js], O, [g], Isl, for 
easy reference and completeness, a brief account is presented here. The INM model is based 
on quantum-mechanical infinite nuclear matter rather than the classical liquid-drop used in 
the traditional Bethe-Weizsacker (BW) type mass formulas. In this model, the ground-state 
energy E^{A, Z) of a nucleus with neutron number A^, proton number Z, mass number A 
and asymmetry /3 = (A^ — Z)/ {N + Z) is considered equivalent to the energy of a perfect 
sphere made up of infinite nuclear matter at ground-state plus the residual characteristic 
energy called the local energy r/. rj mainly consists of energies due to shell, deformation and 
diffuseness etc., which are intrinsic characteristic properties of a given nucleus and can be 
considered as its finger-print. Thus a nucleus possess two categories of properties, namely, 
the global one represented by the INM sphere and the individualistic one by rj (A,Z). So 

E^{A,Z) = E'{A,Z) + r^{A,Z) (1) 

with 

E'{A,Z)=E{A,Z) + f{A,Z), (2) 

where 



f{A,Z) = afA2/3 + a[,[Z2-5(— )2/3z4/3]A-V3 

IGtt 

+ai,A^/^P^ + ai,A^/' - Z) (3) 

denotes the finite size effects and E{A, Z) is the energy of the infinite nuclear matter con- 
tained in the sphere. This sphere is hereafter referred to as the INM sphere, and the super- 
script I stands for the INM nature of the coefficients. Here , a^, a^^ and a^^ are the surface. 
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Coulomb, surface-symmetry and curvature coefficients and S{A, Z) is the usual pairing term, 
given by 

S{A, Z) = +AA^^/^ for even — even nuclei 
= for odd — A nuclei 

= -AA'^^^ for odd - odd nuclei (4) 

Eq. (P ) now becomes 

E^{A, Z) = E{A, Z) + /(A, Z) + r/(A, Z) (5) 

Thus the energy of a finite nucleus is written as the sum of three distinct parts: An 
infinite part E , a finite part / and a local part r] . All these three parts are considered 
distinct in the sense that each of them refers to a different characteristic of the nucleus and 
as such, is more or less independent of each other. Eq. (jS] ) is our required mass formula 
which provides a direct link between finite nuclei to nuclear matter. Its three functions E, f 
and 7] are to be determined. 

The term E{A, Z) being the property of nuclear matter at the ground state, must satisfy 
the generalized [4] HVH theorem given by 

I + P^^^ = ^[(1 + /^)^" + (1 - ^^'P^- (6) 
Using Eq. (jS] ), the INM Fermi energies e„ and ep can be expressed in terms of their 
counterparts of finite nuclei denoted through suffix 'F' as 



e^-(df/dN),-{dv/dN), 



= 



(df/dZ)^ - (dv/dZ)^ (7) 



where = {dE^/dN)z and = {dE^ /dZ)^ • Using Eqs. ([5]) and ([7]), the relation P) 

can be recast as 



^ + ^ = ^[(1 + PVn + (1 - + S{A, Z) (8) 
+\[{^ + P){^)z + {1-P){%)nI (9) 
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where, 

S{A, Z) = f/A - {N/A){df/dN), - {Z/A){df/dZ)^ (10) 

is a function of all the finite-size coefficients , a^, afg and a^^ which are global in nature. As 
noted earlier, the local energy t] refers to a the individualistic characteristic of the nucleus, 
while the global function S{A, Z) refers to the bulk properties which are global in nature. 
Also the above Eq. (|9]) does not have any coupled term involving these two functions i] and 
S. Therefore, we make the ansatz of splitting the above equation into the following two 
equations. 



[(i + /3K + (i-/3)<] + 5(A^) (11) 



or equivalently 



S{A, Z) = ^-1[{1+ (3)e^ + (1 - (12) 



and 



^^^^ = ^[(1 + mdv/dN)z + (1 - mdv/dZU]. (13) 

If individually each of these Eqs. (|T2l [T3l) is satisfied, then the original Eq. (jO]) is satisfied. We 
would like to stress here that the validity of these two equations have been well demonstrated 
a'posteriori in numerical calculations elsewhere js, 5, 6|. 

These three equations define the INM model completely; Eq. (6) determines the finite- 
size coefficients and etc. of the INM sphere, Eq. (7) determines the global parameters 

and while Eq. (8) exclusively determines the local energy t]. Thus once the three 
functions E{A, Z), f{A, Z) and ri{A, Z) are determined, the energy of the nucleus (A,Z) is 
obtained using the mass formula Eq. (3). 



III. CALCULATION 



It should be emphasized here that the above three Eqs. (6-8) exclusively determine the 
surface, the bulk and the local energies of a given nucleus separately. The decoupling of these 
three quantities provides the necessary breakthrough 8|, |9| for the correct determination of 
their respective parameters inhibiting cross correlations amongst them. This ledjs, 9| to the 
resolution of the tq — paradox and determination of the saturation properties including the 



incompressibility of nuclear matter entirely from nuclear masses. The parameters a^, a^, ajg 
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and a^^ characterizing the INM sphere are global in nature. The surface symmetry a^^ term 
cancels to the extent of 66% due to particular combination of data occurring in it. Further 
being a higher order term, it is quite weak compared to the principal term. Combination of 
these two features makes it redundant and rather symbolic at numerical level in the INM 
model thereby playing no significant role. We have demonstrated this through exhaustive 
numerical calculation in references [js, 9|]. The same is true for the curvature a^^ term. So we 
have dropped these two terms. We would like to emphasize that in the INM model we have 
adopted the principle to determine all the parameters using exclusively the ground-state 
masses of nuclei only as the model is intended to predict the ground-state masses. In the 
BW-like mass models 16|,ll7|] the the value of surface-symmetry coefficient agg, is usually fixed 
using properties other than the ground-state masses like fission barrier etc. and the curvature 
term is normally dropped as its coefficient cannot be determined uniquely. It is true that by 
dropping the surface-symmetry term in our model, isospin-dependence is exclusively taken 
care of by the bulk symmetry term only and any left-over isospin dependence arising from 
the surface-symmetry of a given nucleus is taken care of by the characteristic local energy 
ri{N, Z). Thus this in no way denies any loss of generality regarding explicit absence of the 
isospin-dependent surface-symmetry term in the INM model. 

We have also made an exhaustive study to include other higher order terms like Nolen- 
Schiffer charge- asymmetry and proton-form factor etc., however they could not be deter- 
mined uniquely using the masses and the neutron- and proton-separation energies in the 
framework of our model as they exactly cancel out in the INM Eq. 6 (see Ref.j^, 9| for 
details). In BW-like mass models the coefficients of these higher order terms are fixed using 
properties other than the ground-state masses. They vary widely from one model to another. 
For example, in the Finite- Range Droplet (FRDM) modeljl^, the value of charge-asymmetry 
coefficient is 0.436 MeV while in the Finite- Range Liquid-Drop (FRLDM) model 16| its value 
is 0.10289 MeV. In view of such wide variation in the value of this parameter, we feel in 
general, that it is not desirable to assign some values to them from external considerations 
using properties other than the masses. Thus although their explicit representation in the 
INM model is untenable, however their collective effect is implicitly accounted for through 
the local energy eta being determined by solving Eq. (8) using known experimental masses. 
Thus we have altogether four global parameters a^, a^, a^, al and the pairing coefficient 
A. Since these parameters are universal in nature being valid for all nuclei, they have been 
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meticulously determined earlier {g, |8|, |9| taking experimental masses with error bar less than 
60 keV to avoid pollution from other available data with higher error. With these criteria, 
the masses of nuclei in the valley of stability were chosen together with their neutron and 
proton separation energies in fitting Eq. (6) and Eq. (7). Although in the present mass 
table we are using extra 314 nuclei for mass predictions compared to our 1999 mass table, 
they are not necessarily suitable for the determination of the global parameters as they are 
more exotic in nature. So we have adopted the same values for the global parameters as 
determined for the 1999 mass tablejojls, 9|. These values are 

al = 16.108 MeV 



a, 



A 



24.06 MeV 
19.27 MeV 
0.7593 MeV 
11.505 MeV. 



Once the global parameters defining E{A, Z) and f{A, Z) are known, the empirical values 
of 7] for all known nuclei are determined using the mass formula Eq. ([5]) . Prediction of 
masses of nuclei in the unknown regions amounts to predicting the ?7's of such nuclei, rj 
satisfies Eq. (fT3|) which is a partial differential equation in (N,Z) space. Being expressed in 
terms of rate of change of t] with N and Z, it is endowed with good predictive ability like well- 
known equations in Physics. With proper initial conditions given by experimental values 
of 7] in known regions, long range prediction of 77s can be obtained in unknown regions by 
solving the Eq. (fT3ll . In all our previous publications on the INM model we have presented 
in detail the construction of the grid in (N, Z) space and the method of solution, which we do 
not feel any necessity to repeat here. However, we would like to reiterate only the ensemble 
averaging process which was introduced for the 1999 mass table predictions 6|], which in fact 
was crucial for getting proper solution of the rj-Eq. (8) . Fixing a hexagonal grid for 58 nuclei 
in a given region in (N,Z) space (see Fig. 1), Eq. (8) is solved to obtain definite values of 1] 
for those nuclei. Shifting the grid in all possible directions in (N,Z,A) space, an ensemble of 
about 70 alternate values for a given nucleus is generated. The most probable value is then 
obtained by performing ensemble-averaging using the standard Gaussian weighted method 
given by 
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FIG. 1: The vertices in the hexagon designated by neutron and proton numbers show possible 
connections among the nuclei taking part in the solutions of the local energy //-equation. (jl3p . 

^ ^ EiVi exp[-{{ai - ao)/armsf] ^^^^ 
exp[-{{ai - ao)/armsf] 
where (Tq is the least of all the CTjS for a given nucleus and arms is the rms deviation of all 

the 2198 known nuclei used in the fit. Finally we use the global terms E{A, Z) and f{A, Z) 
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in Eq. (3) along with the value of rj so determined to predict the mass of a given nucleus in 
any part of the nuclear chart. 




Neutron Number N 



FIG. 2: Differences between the experimental (Ref.[l8]) and calculated binding energy(BE) for 
2198 known nuclei. 



As before we observe a systematic deviation in the binding energies of N=Z nuclei as a 
function of mass number as shown in the Fig. 2 of reference {gI . There are altogether 31 
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such nuclei. The binding energy systematics resemble very closely the Wigner contribution 
usually taken into account in the droplet-like models. As before we have taken the Wigner- 
like term given by W = Wq/A5n,z as a correction to be applied for such nuclei. We have 
taken the value of Wq as 32 MeV as determined before 6|. 

Once the global parameters and the local energies are obtained, the INM mass Eq. (3) 
is used to calculate the nuclear masses for the entire nuclear chart. Thus the present mass 
table is constructed comprising 6727 nuclei in the range 4 < Z < 120 and 8 < A < 303 
extending up to the drip lines. It may be noted here that our present mass predictions are 
confined up to the drip lines and in most cases three to four steps beyond the drip lines 
in the nuclear landscape. We feel that this coverage of the nuclear landscape is sufficient 
enough for most applications both in Nuclear Physics and Astro-Physics. That is why the 
present mass table has 6727 mass predictions compared to 7208 in our 1999 mass table 6|. 



The rms deviation for the fit to 2198 data-nuclei 
deviation is 1.3keV. 



18| is found to be 342 keV while the mean 



IV. RESULTS AND DISCUSSIONS 
A. General Features 

To check if there is any local fiuctuations in the quality of our fit, we have plotted in 
Fig.2 the differences between the experiment al[18i] binding energies (BE) and the fitted ones 
termed as BE residuals, in the different regions of N — Z plane. The various symbols in the 
figure represent the magnitude of the BE residuals. It is satisfying to see that the magnitude 
of the residuals mostly lie within 350 keV in accord with the rms deviation of 342 keV. In 
Fig. 3, we have also plotted the BE residuals as a function of mass number A to study if 
any residual systematics are present, which would signal the deficiencies of our model. It is 
interesting to see that the residuals are evenly distributed around the A-axis mostly lying 
within 350 keV, which is well corroborated by the small mean deviation of 1.3 keV. 

We would like to point out below the unique features of this mass formula which endows 
it with such low rms and mean deviation. 

(i) The INM model uses not only the nuclear masses like other mass models, but also 
additional two sets of data namely the energy differentials {dE/dN)z and {dE/dZ)^. 
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FIG. 3: The deviation as a function of mass number A of the calculated binding energies (BE) 
from the experimental data taken from Ref. [l^for 2198 nuclei. The two horizontal lines parallel 
to A-axis drawn at ±0.35 MeV enclose most of the BE residuals. 

It may be mentioned that for any function in general, all the derivatives at a point 
are its independent properties apart from the value of the function at that point. 
Thus effectively the model uses three times the mass data pertaining to three different 
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properties of nuclei unlike other models for fitting its parameters. 

(ii) The model uses the same set of data two times separately in solving Eqs. (6) and (7). 

More over in solving Eq. (8) it uses several tens of times the same set of data. 

(iii) Unlike the droplet-like mass formulas, the global parameters al and ajj are determined 
by one equation (Eq. 7) and the and etc. are determined by another equation 
(Eq. 6), which allows no cross-correlation among them thereby leading to correct 
values of these parameters. 

(iV) Since HVH theorem is valid for three-body and other multi-body forces, the INM 
model being based on it takes them implicitly into account. This is very important 
since three-body forces have been recognized to contribute significantly to nuclear 
structure and the saturation properties. 

In our least-square fit, we have not included the masses of those nuclei given in the mass 

n 

Table of Audi and Wapstra[l8] which are not experimentally measured but estimated from 
systematic trends (referred to as Audi-Systematic Data). These nuclei mostly lie in the 
vicinity of known mass surface with relatively large error bars. There are altogether 936 
such nuclei. It may be of interest to compare our predictions with these data as presented in 
Fig. 4. It is satisfying to note that our predictions compare well with the Audi- Systematic 
Data. 

B. Predictive Potential of the INM model 

We would like to critically examine the predictive potential of the INM model basing on 
two aspects; namely the structural nature of the local energy t] -systematics, and comparative 
global analysis of the predictions of different mass models. 

1. Structural Nature of r]- systematics: 

The key element in the INM model is the local energy rj which embodies all the character- 
istic properties of a given nucleus, namely the shell effect, the deformation, the diffuseness 
etc., and possibly other unknown properties. Therefore it is expected that study of rj should 
reveal some important features of nuclear dynamics. In such a study[l^ the plot of ri{A,Z) 
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20 40 60 80 100 120 140 

Neutron Number N 



FIG. 4: Differences between binding energies (BE) of the Audi-Systematics data (Ref. la]) and 
the calculated values of 936 nuclei. 

as a function of N for a given Z shows Gaussian peak at shell-closures in the valley of stabil- 
ity. Such plots using the present empirical values of rj for all the known nuclei are presented 
in Fig. 5, where well-defined Gaussian structures at the magic numbers 50, 82 and 126 
are prominently seen. This feature is suggestive that the r^-systematics in the unknown 
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regions may reveal new shell-closures in the neutron-rich region close to drip-line which is 
corroborated in the following study. 




FIG. 5: Empirical values of rj as functions of N for all Z showing structural nature of r/-isolines. 
Vertical lines drawn at the Magic Numbers 20, 50, 82 and 126 and the semi-magic number 28 
indicate Gaussian structure. 



In 2004, we ll| have shown the predictive potential of the INM model by predicting new 
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neutron magic numbers 100, 150,164 and proton magic number 78 and weakly proton-magic 
numbers 62 and 90 along with new islands of stability around ^^^Sm,^^® Pt and "^^^Th in the 
drip-line regions. This prediction signals a new phenomenon where shell effect overcomes the 
instability due to repulsive components (triplet-triplet, singlet-singlet) of nucleon-nucleon 
interaction, which is complementary to the phenomena of super heavy elements and fission 
isomers, in which repulsive coulomb instability is overcome by the same. While the latter 
elongates the stability peninsula the former broadens it. This prediction was done on the 
basis of and rj systematics obtained in our INM mass predictions. These results were 
corroborated through microscopic studies in RMF theory and Strutinsky shell correction 
calculations. Here we reinvestigated the same with our new mass predictions. In Fig. 6, we 
have plotted S2n in the upper panel, rj in the middle panel and the S2n- differential given 
by S2n{N, Z) — S2n{N + 2, Z) in the lower panel, which reaffirms our earlier result stated 
above. Thus the ability of the INM model for long range predictions looks promising. 



2. Comparative Global Analysis of Mass Models: 



Presently about 2200 nuclei are known and masses of another 5000 to 7000 nuclei have 
been predicted by different mass formulas. Some of them are expected to be synthesized in 
the RIB facilities in the comin g ye ars. Therefore the predictive ability of different mass for- 
mulas is under serious scrutiny 12Nl4| in the recent years. For this purpose, the predictions 



of different mass models 



19| for Sn isotopes with neutron number ranging from 45 to 110 



have been compared with those of Duflo-Zucker 
ously low rms value of 375 keV. Such studies [12 



as the reference chosen for its conspicu- 



14| reveal that all mass models 16|, |l9| show 



good agreement in the known region in the valley of stability, however the predictions di- 
verge randomly without showing any correlation as one moves away to the unknown regions 
on either side where experimental results are not available. This has raised serious question 



about the efficacy of mass formulas 
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19| as the tools for prediction. For our discussion 



here we carry out such comparisons more extensively with respect to Duflo-Zucker, RMF 
and our present INM model predictions not only for Sn isotopes, but also for Pb and Ca 
isotopes. Here for a nucleus (N,Z), the BE difference BEDi{N, Z) in the model i is obtained 
as 

BEDi{N, Z) = BE,{N, Z) - BEref{N, Z) (15) 
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FIG. 6: S2n,'n and S'2n-differential (see text ) showing occurrence of new shell-closures in 
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where BEi{N, Z) and BErefiN, Z) are the predicted BEs in the mass model i and the 
reference mass model respectively. With Duflo-Zucker taken as the reference model, results 
are plotted in Fig. 7 for the Sn, Pb and Ca isotopes separately. It can be seen that there 
is unanimity of all the mass models {igI, 19 | on good agreement with experiment in the 
known region in the valley of stability; however the predictions diverge as one moves away 
to unknown regions on either side. Since the masses of the known regions have been used 



in the fit by all the mass models 
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19|, the agreement with data is to be expected. What 



is worrying is in the unknown regions where they do not exhibit a common trend like rising 
and falling, but diverge without any correlation. The divergence with the same intensity is 
also seen when the model of Moeller and Nix et al[16| is used as the reference model in Eq. 



about the efficacy of mass models 



(10) as is seen in Fig. 4 of Ref. 12 1. This uncorrelated divergence has raised the question 
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19| as a whole 12 -[ij]. One is constrained to think 



whether development of different mass models is a parameter game only without potential 
for reliable predictions! A probable explanation of this phenomenon may be as follows. 

It is a common feature with most mass models that the degree of success varies from 
region to region in (N,Z) space even in known domains and may be more so in unknown 
domains. It is quite likely, the degree of accuracy of the predictions of a given model may 
not match with the corresponding ones of the reference model in the same (N,Z) domain 
resulting in randomness with no common trend, which appear as divergences seen in Fig. 7. 

A microscopic study based on nuclear Hamiltonian may throw light on the above issue. 
We feel masses calculated in Relativistic Mean Field (RMF) theory can qualify to serve as a 
reference mass model. The popular parameter set widely used in literature with conspicuous 
success is the NL3 set 2l| . The 7 parameters of this Lagrangian have been determined by 
fitting the data of ten nuclei only, and with the interaction so determined the predictions [22 1 
of 1315 even-even nuclei give rms deviation of 2.6 MeV. It must be noted that the masses of 
these nuclei were not fitted unlike other models. So they should be considered as predictions 
in the unknown region. Therefore similar validity of this interaction in the framework 
of the RMF for a few hundred nuclei in the true unknown regions may be expected in 
a somewhat lesser scale. The BEDi's calculated with the RMF predictions as reference 
for all the isotopes of Ca, Sn and Pb are shown in Fig. 8. It can be seen that for each 
chain of these isotopes, the predictions in the unknown regions show a common trend in 



general for all the mass models 16 



19| . although qualitatively they differ from one another. 
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FIG. 7: Binding energy difference BEDis{ Eq. 10) in Ca,Sn and Pb isotopes of various Mass 



Formula[l6l. Il9l | predictions with respect to those of Duflo-Zuker 



201 1 as the reference. 



Therefore we may conclude that the predictions of mass models in the unknown regions are 
not necessarily arbitrary or random but show definite trends, which therefore can be used as 
useful guidelines for theoretical and experimental studies. So RMF mass predictions provide 
a uniform substratum for the whole (N,Z) domain suitable to serve as a reference. 
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FIG. 8: Binding energy difference BEDis{ Eq. 10) in Ca,Sn and Pb isotopes of various Mass 
Formula [l^. \\\ predictions witfi respect to those of RMF[2_2:] as tlie reference. 

Now in the context of our present mass predictions in the INM model, we use it as the 
reference model in Eq. (10) to calculate the BEDiS for all those isotopes of Ca, Sn and 
Pb. The results so obtained are presented in Fig. 9. It is satisfying to find that no random 
divergences appear in the unknown regions, rather the expected common trends are seen in 
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all the three cases. Thus like the RMF mass predictions, INM mass predictions can also be 
used as substratum for suitable reference purposes. This implicitly bears out our expectation 
about the predictive potential of the INM mass model. 
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FIG. 9: Binding energy difference BEDis{ Eq. 10) in Ca,Sn and Pb isotopes of various Mass 
Formula [ifil. \\\ predictions with respect to those of the present INM predictions as the reference. 
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C. Shell Quenching in = 82 and N = 126 Shells 
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FIG. 10: Shell gaps given by S'„(Ar^ -1,Z) - SniNm + 1,2) for iV^=82 and 126 for various 
proton number Z showing shell-quenching in N=82 and 126 shells. 



The shell structure is a very fundamental feature in nuclei. It plays a crucial role in 
building of different models for their description. In the valley of stability occurrence of 
well-defined shell structure in the form of magic numbers is well established and extensively 
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studied. In the INM model the same feature is also seen in the systematics of local energy 
1] when these are plotted as isolines as shown in Fig. 5. One can clearly see such sharp- 
peaked Gaussian structures at the neutron magic numbers 50, 82 and 126. In general, 
whether such features persist in the unknown regions far from stability is of great current 
interest. Mass measurements of neutron-rich nuclei in the light and medium heavy regions 
have shown 23(] that = 20 shell gap vanishes for nuclei with Z < 12. Such shell quenching 



has also been seen for the = 50 shell 2M- hoped that this phenomenon should also 
be seen for = 82 and 126 shells. However conclusive proofs presently are lacking as 
masses of the relevant nuclei lying in the drip-line region have not yet been measured in the 



laboratory. 



Analysis of the abundances of various elements in r-process nucleosynt 



hesis is 



Id, 
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suggestive [251] of the quenching of A^ = 82 shell. Our 1999 mass predictions show 
quenching in A^ = 82, and A^ = 126 shells. The present mass predictions reaffirm the same 
more convincingly. In Fig. (10), the shell gaps calculated using the present mass predictions 
are shown for A^ = 82 and A^ = 126 shells. Gradual vanishing of the gaps in both the cases 
are seen as one moves towards the drip-line regions. 

D. Islands of Inversion 



Neutron-rich nuclei lying in the proximity of the n-drip line are expected to show some 
exotic features in nuclear structure with dramatic departures from the normal trends ob- 
served in the /3-stable valley. A well-known phenomenon having such a feature termed as 
'island of inversion' has been observed [22] in neutron-rich nuclei around N=20 region quite 
early, which manifests in enhanced binding of those nuclei centering around ^^Na. It has 
been always a challenge to mass formulas and other structure models to explain this phe- 



nomenon. Extensive theoretical and experimental studies 28|, 1291] carried out over the years 
have concluded that the N=20 shell-closure in this region is broken by the intruder states 
from the p/-shell, inducing strong collectivity giving rise to high deformation. This has 
resulted in enhanced binding for Na and nearby nuclides in this region. Therefore with the 
present mass predictions in the INM model, we have attempted to make a through search 
throughout the nuclear landscape to see if our masses can account for the above known 
island of inversion and also identify possible new such islands. With this view we have 
plotted the S2n isolines for all the nuclides from Z=7 to 105 in Figs. 11-13. The typical 
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FIG. 11: The two-neutron separation energies S2n as function of N for nuclei with Z= 7 to 36 
drawn as sohd, broken, dash-dotted etc. hnes. The corresponding Experimental [3] values where 
ever available are indicated by various symbols to distinguish the neighbouring isolines. 

sharp fall of at the shell-closures clearly reproduce the well-known magic numbers 8, 
20, 50, 82 and 126 in conformity with experiment. However the monotonic decrease with 
increase of neutron numbers N in the /3-stable valley gets arrested in the neutron-rich region 
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for Z=10 and 11, agreeing with the observed island of inversion around '^^Na. Thus our 
mass predictions reproduce this feature reasonably well. 




FIG. 12: Same as Fig. 11 but for Z= 36 to 70. 

It can also be seen in Fig. 11 that there is a region spanned by Z=17 to 23 and N=38 to 42, 
where isolines exhibit the same feature of enhanced stability suggesting the existence of 
another island of inversion. This may originate due to breaking of N=40 shell by the intruder 
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states from the sdg shell and thereby inducing strong deformation. It is indeed satisfying 
to note that in the recent years strong evidences both theoretical and experimental have 



emerged 
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3l| supporting the existence of this island of inversion centering around ^^Tz. 



These two islands of inversion do suggest that these two may not be the isolated cases, and 
this phenomenon may be a general feature of nuclear dynamics in the exotic neutron-rich 
regions close to n-drip line, where breaking of shell- closures by intruder states from higher 
shells are highly plausible. In fact our extensive systematics in the high-mass region 
presented in Fig. 12 reveal two more islands in the heavy-mass region delineated by Z=37- 
40, N= 70-74; and Z=60-64, N= 110-116, where these may be due to breaking of N=70 and 
N=112 shells by the intruder states from the pfh and sdgi shells respectively. 
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FIG. 13: Same as Fig.ll but for Z= 70 to 105. 
V. CONCLUSION 

The Infinite Nuclear Matter model of atomic nuclei developed to its full potential over the 
last thirty years has been presented. The model has uncovered a new facet of the nucleus in 
which the nuclear matter characteristic has been revealed to be predominant. The abiding 
picture since the emergence of nuclear physics in 1930s is that the nucleus can be considered 
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as a classical liquid drop which formed the foundation of Bethe-Weizsaecker(BW) mass 
formula and its more refined version like the droplet model etc. These liquid drop based 
models are traditionalally used to define infinite nuclear matter and its saturation properties, 
determinable from nuclear masses. This has lead to many ambiguities and inconsistencies 
like ro-paradox, indeterminability of incompressibility etc. In contrast, the INM model has 
replaced the classical liquid by a quantum mechanical many-fermionic liquid which is the true 
nature of nuclear matter composed of the nucleus, by founding it on the generalized HVH 
theorem of many-body theory. Thus the liquid and the shell features, the two main features 
of nuclear dynamics could be taken into account non-perturbatively at the microscopic 
level. The model has succeeded in scpcrating the energy of the nucleus into a universal 
part represented by a sphere made up of infinite nuclear matter common to all nuclei, and 
a charecteristic part r} called the local energy which varies from nucleus to nucleus. The 
universal part is described by only four parameters charecteristsic of the nuclear matter and 
77 being determined from nuclear masses through the solution of a differential equation of 
the model. 

The effect of three-body has been included inadvertently by the use of HVH theorem, 
which is valid for such forces. Because of this proper physics input, the saturation density 
of nuclear matter could be obtained from nuclear masses in agreement with the value de- 
termined from electron-scattering on heavy nuclei, and thereby resolving the long standing 
ro-paradox. The energy per nucleon E/Aoi nuclear matter at ground-state has been con- 
sistently obtained from nuclear masses also. Thus both the saturation properties of nuclear 
matter E/A and poo are obtained from the same source using the INM model. The third 
property of INM, namely the Incompressibility could be determined in more or less model 
independent manner using the same mass data that includes nuclear masses, and neutron 
and proton Fermi energies. Thus the data used for the determination of all the three sat- 
uration properties numbering about 4000, are well- respected as the most precision data in 
nuclear physics and that too pertaining to ground-state only. 

The mass formula based on the INM model has the unique feature that it is formulated in 
terms of differential equations and thereby imbued with the ability for prediction of masses 
in unknown region. The root mean square deviation and the mean deviation of the mass-fit 
to 2198 known masses yields respectively 342 keV and 1.3 keV only the lowest in hterature. 

The model has introduced in nuclear physics a new entityi, namely, the local energy 77. It 
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contains all the characteristic properties like shell, deformation, diffuseness etc and therefore 
exceeds in scope compared to the usual shell-effect introduced by Strutinsky. In known 
regions, it is determined from experimental mass, and in unknown regions it is predicted by 
solution of a differential equation using the known ones as boundary conditions. The former 
being unique is endowed with the predictive potential, whereas the latter depending on the 
chosen meanfield is deficient in this respect. The systematics of r] reveals new magic numbers 
in exotic neutron-rich nuclei in the dripline regions of the nuclear chart. This predictive 
ability is not contained in the usual shell effects, the latter being dependent on mean field 
which varies from interaction to interaction. The systematics of r) predicts three islands 
of stability with N=100, Z=62; N=150, Z=50 and N=164, Z=64 around ^^^SM, ^^^Pt 
and '^^'^Th supported through microscopic studies in the framework of RMF theory and 
Strutinsky shell-correction calculation. This reveals a new effect in nuclear physics, where 
shell overcomes the instability due to repulsive components (triplet-triplet, singlet -singlet) 
of nuclear force, analogues to the superheavy elements where repulsive Coulomb instability 
is stabilized by the same, while the latter effect elongates the stabihty peninsula the former 
broadens it. 

In comparative global analysis of all the mass models, the unique predictive ability of 
the INM mass model has been established by using its predicted binding energies as the 
reference and studying the systematics of the predicted masses of other mass models. 

The shell-quenching in N=82 and 126 shells have been found using the INM masses, much 
anticipated but not predicted by any other mass models 

The model has also shown the existence of new islands of inversion dehneated by Z=37- 
40, N=70-74; Z=60-64, N=110-116; thus showing that the phenomenon is quite general and 
not confined to the well-known islands of inversion around ^^Na for a long time in nuclear 
physics. Thus in this respect, the INM model has been singulalarly successful. 

The unique features of the INM mass model are: 

i) It does not use any effective interaction in any way, and thereby it is free from the 

vagaries of the uncertainty and from region to region variation. 

ii) It uses only the nuclear masses and Fermi energies determined solely by experiment. 

iii) It takes into account three-body, and other multibody forces (if present) in nuclei. 
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4) It is built in terms a differential equation and thus equiped with physics for prediction. 
The distinctive success of the INM model may be attribited to these features. 
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VI. EXPLANATION OF TABLE 



VII. TABLE 



Calculated Nuclear Ground-State Binding Energies and Mass Excesses, Com- 
pared to Experimental Mass Excesses Where Available along with One-neutron, 
Two-neutron, One-proton and Two-proton, a-particle Separation Energies 

Z Proton number. The mass table is ordered by increasing proton number. 

The corresponding name of each element is given in parenthesis. 

Neutron number. 
A Mass number. 

BE(MeV) Calculated binding energy B(N,Z) of a nucleus (N,Z) using Eq. (3). 
ME(MeV) Calculated mass excess. Nuclei that are unstable against 

one-nucleon and two-nucleon decay are denoted by symbols (f) and (J), respectively. 
MEexp(MeV) Experimental mass excess of Audi-Wapstra[18]. 
Err (Mev) The error associated with the Experimental mass excess[18]. 
SniMeV^ Calculated One-neutron Separation Energies given by B(N,Z)-B(N-1,Z). 
S'2n(MeV^) Calculated Two-neutron Separation Energies given by B(N,Z)-B(N-2,Z). 

(MeV) Calculated One-proton Separation Energies given by B(N,Z)-B(N,Z-1). 

(MeV) Calculated Two-proton Separation Energies given by B(N,Z)-B(N,Z-2). 
5a (MeV) Calculated a-Particle Separation Energies given by B(N,Z)-B(N-2,Z-2)-28.296. 



